A dynamical study of fission process and estimation of prescission neutron multiplicity 
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The dynamics of fission has been studied by solving Euler-Lagrange equations with dissipation 
generated through one and two body nuclear friction. The average kinetic energies of the fission 
fragments, prescission neutron multiplicities and the mean energies of the prescission neutrons have 
been calculated and compared with experimental values and they agree quite well. A single value 
of friction coefficient has been used to reproduce the experimental data for both symmetric and 
asymmetric splitting of the fissioning systems over a wide range of masses and excitation energies. 
It has been observed that a stronger friction is required in the saddle to scission region as compared 
to that in the ground state to saddle region. 



I. INTRODUCTION 



The studies of fission dynamics has become a subject of current interest due to availability of recent experimental 
' data for prescission neutron multiplicities in heavy ion induced fusion-fission reactions [l|-|lC|]. From the analysis 
of the data, it is by and large established that standard statistical theory could not account for the observed large 
multiplicity of the prescission neutrons. This discrepency is, nowadays, believed to be of dynamical origin, and thus 
led to enormous activities to understand theoretically the dynamics of fission process j|TT| |22| . 

Dynamics of fission consists in the study of gradual change of shape of a fissioning compound nucleus. The shape 
is globally characterised in terms of elongation parameter, the neck radius and mass ratio of two fragments. These 
variables are usually referred to as collective variables and dynamics is understood in terms of the evolution of these 
collective variables. It has been observed that time taken by the compound nucleus from its state of formation to 
t~^. . the state when it gets dissociated into the two fragments is 'too short' when the shape evolves under the effect of the 
0^ 1 conservative forces derived from the Coulomb and surface energy of the compound nucleus. With the average neutron 
evaporation rate determined by excitation of the compound nucleus, such large prescission neutron multiplicity data 
■ can not be explained in any manner with this 'short' time scale of fission process. This feature leads one to believe 
that the dynamical process of shape evolution gets inhibited for a considerable amount of time. This retardation in 
the dynamics of collective variables is effected through a mechanism of dissipation by assuming the fissioning nucleus 
as a liquid drop. The dissipation is usually realised by introducing a friction term in the dynamics of shape evolution. 
Thus, the explanation of the observed large neutron multiplicity points to two important features of the fission process; 
the time scale of fission and the viscosity of the nuclear fluid. The dynamics of fission is then picturised as a dissipative 
process where initial energy of the collective variables get dissipated into the internal degrees of freedom of nuclear 
fluid giving rise to the increase in internal excitation energy, which, in turn, is responsible for the evaporation of 
' prescission neutrons. 

There have been different approaches to study this problem. One approach is to solve the Langevin equations for 
the collective variables flll|Jl2|]l^ , |l5|| . In this approach, one assumes the collective variable as the 'Brownian particle' 
interacting stochastically with large number of internal degrees of freedom constituting the surrounding 'bath'. The 
systematic dissipative force is assumed to be derived from the random force averaged over a time larger than the 
collisional time scale between collective and internal degrees of freedom. The random part is modelled as a gaussian 
white noise which causes the fluctuation of the physical observables of the fission process such as kinetic energies, 
yields of the fission fragments etc. In the other approach, one solves the multi-dimensional Fokker-Planck equation 
]Tg| Q, which is a differential version of Langevin equation. In both of these approaches, the calculated average 
kinetic energies, average yield of the fission fragments and neutron multiplicity compare more or less well with the 
respective experimental data. 

Following an alternative approach, we proposed a dynamical model of fission |^,^| , which could explain fairly well 
various features observed in fusion-fission reactions of lighter systems. In that model, we solved the Euler-Lagrange 
equations for the collective variables instead of solving the Fokker-Planck or Langevin equations. The fluctuations are 
introduced at the initial level of the dynamics by random partitioning of available energy of the 'nascent' compound 
nucleus between the collective and the intrinsic degrees of freedom and attributing the former to the generator of 
the displacement of the collective variable as an initial condition of Euler-Lagrange equation. The random initial 
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momentum given to the collective variable gives rise to different trajectories. Some of these trajectories would cross 
the fission barrier and have the 'fission fate'. The main difference between our approach and those of the others is 
that in our approach the dynamics is deterministic with randomisation at initial level, while in other approaches, the 
dynamics is mainly stochastic in nature. However, in all these approaches the dissipation is introduced in terms of 
frictional forces. 

In our earlier calculations |2l|]22[ | we assumed a schematic shape of the fissioning system which comprised of two 
leptodermous spheres connected by a cylindrical neck. This particular idealisation was introduced by Swiatecki [23]] 
to simplify the calculation of various ingredients such as conservative forces etc. After the saddle point, the neck 
gets constricted and at the scission point the spheres get detached forming two fission fragments. The calculation 
with such simple shape produced results which agree quite well with the experimental data, particularly for the nuclei 
lying below the Bussinaro-Gallone point. It is, however, known that for heavy nuclei lying above the Businaro-Gallonc 
point, the fission shapes are highly deformed, and elongation is quite large before the nucleus reaches the scission 
point. Obviously, the simple schematic shape parametrisation [^3j used earlier would not be expected to mimic such 
large deformed shapes. Besides, the dissipation of collective energy to the internal nucleonic degrees of freedom could 
happen through two different mechanisms. One is due to two-body collision of the nucleons inside the nuclear fluid, 
the other is due to collision of nucleons with changing surface of the nuclear fluid, which is more commonly known as 
one-body or 'Wall' friction and both of them depend very sensitively on the surface profile of the fissioning system. 
Hence, with these considerations in mind, in the present paper, we have used a more realistic shape parametrisation of 
the heavy fissioning systems to study the temporal evolution of fission shapes, and, subsequently, to calculate various 
observables of the fission process, such as prescission neutron multiplicities, total kinetic energies etc. 

The present paper has been arranged as follows. In Sec. || we describe the model and the statistics used in 
calculating the aver age kinetic energies, prescission neutron multiplicities etc. The results of the calculation are 
discussed in Sec. III. Finally, concluding remarks are given in Sec. |TV|. 



II. THE MODEL 

The present model is a generalisation of the schematic model developed by us [^TJ to study the fusion-fission process 
for light nuclear systems lying below the Businaro-Gallone point. The schematic parametrisation for the evolution 
of fission shapes used in the earlier work has been replaced by a more realistic parametrisation for the same that is 
applicable for heavier nuclei with fissility parameters above the Businaro-Gallonc value. 

A. The shape 

The shape of the nuclear surface is assumed to be of the form 

p 2 {z) = c~ 2 {c 2 - z 2 )(A + Bz 2 + azc), (la) 
where the coefficients A and B are defined as 

A = c~ 1 -B/5, (lb) 



B=(c-l)/2. (lc) 

This is a specific form of the surface introduced by Brack et. al. [Q. The quantity c corresponds to the elongation 
and the quantity a is a parameter which depends upon the asymmetry (a sym ) defined below. We may note that the 
surface cuts the z-axis at z — ±c, so that the surface to surface separation along the elongation axis is 2c. There are 
three real points at which the derivative of p vanishes; two of them correspond to maximum of p and between these 
two maxima one minimum occurs at z = z m i„. The portion of volume contained in z — — c to z = z m i„ is defined as 
left lobe and remaining portion of the volume is referred to as right lobe. Asymmetry (a asy ) is then defined as 

_ (A R - A L ) 

a asy — . \ z ) 

A-CN 

where Aqn is the compound nucleus mass, and An, Al correspond to the masses of the right and left lobes, 
respectively. The parameter a is related to the asymmetry a asy through the following relation 
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.11937aL, + .24720<z osy . (3) 



As the shape changes gradually, the coordinates of the two maxima and that of the minimum change. The scission 
point is defined when the minimum point touches the z-axis and it is given by 

c 2 a 2 

•4-^ = 0. (4) 

Therefore, the value of c at which scission occurs depends on a and the dependence is given by 

c sc = -2.0a 2 + . 032a + 2.0917. (5) 



B. The dynamics 

The dynamics is studied by calculating the semi-classical fission trajectories. The trajectories are obtained by 
solving the Euler-Lagrangian equation , 

.. L 2 . d{V c + Vn) (r v 

fir -=- 7r r , 6a) 

pr a or 

hdi = 7* [52 (0a -0) + gi (0i - 0)]g u (6b) 

h92 = lMe 2 -6) + gi(0\ - 9)]g 2 , (6c) 

L = -{i l e l + i 2 e 2 ). (6d) 

The quantities Vc, Vn represent the Coulomb and nuclear interaction potentials and j r) 7 t are the radial and 
tangential components of friction, respectively. I\ , 1% are the moments of inertia of the two lobes and L refers to 
the relative angular momentum. g\ and g 2 are the distances of the centres of mass of the two lobes from the centre 
of mass of the composite dinuclear system and the term [52(^2 — 9) + gi(9i — 9)] represents the relative tangential 
velocity of the two lobes EJ. The variable r is defined as the centre to centre distance between the two lobes. From 
the generalised shape given by Eqn. [la], we first construct the centres of mass of left and right lobes, and call them zi 
and z r respectively. Then r is defined as 

r = \zi-z r \. (7) 
The reduced mass parameter fj,, is obtained from the calculated masses of the two lobes. 

For the non-conservative part of the interaction, we would consider viscous drag arising not only due to two body 
collision but also due to the collisions of the nucleons with the wall or surface of the nucleus. Hence 7r in Eqn. |6a| 
contains two parts; "fJ B and 7 ^? s , for two-body and one-body dissipative mechanisms, respectively. Assuming the 
nucleus as an incompressible viscous fluid, and for nearly irrotational hydrodynamical flow, 7^T B is calculated by use 
of the Werner- Wheeler method [l9|j2C|] and is given by 



where 



-;' i! = »//„ II .,. JUh-iih-) I f dzp 2 [3A' c 2 + \p 2 K 2 ] (8a) 



^(*) = "i^ fdz'p 2 {z>). (8b) 

The quantities A' c , A" are the first and second derivatives of A c (z) with respect to z. po is the two body viscosity 
coefficient. The factor /(ff ) is taken to be 

where x — r/RcNi Rcn being the radius of the compound nucleus. This factor is a consequence of the rotational 



symmetry of the shape ( la) around the elongation axis. The variable r is already defined and the relationship between 
c and x is found to be 
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c = px 2 + qx + f(a), (10a) 
where, p = —.15901, q = 1.03749, and f is given by, 

r = -1.228a 2 - .01896a + .45956. (10b) 
The tangential friction is calculated using the following relation pj]] , 

7 TB = { d<_ fi T B (n) 

The quantity n refers to the neck radius of the composite shape given by Eqn. It is defined as the value of p 
where p 2 has a minima. The variation of n with c for different values of the parameter a is shown in Fig. ^. It is 
evident from the figure that for all values of c, the corresponding values of n are nearly independent of a, and n is 
found to be related to c by the following relation, 

n = -1.047c 3 + 4.297c 2 - 6.309c + 4.073 (12) 

One body dissipative force, -Fdis, is obtained from the rate of energy dissipation, Edi S , by 

d 

Fdis = ~-QrEdis{x) (13) 

where x refers to the rate of change of x with respect to time and Edi S {x) is the rate of energy dissipation at x 
given by 

Edis = 7^p m v <j> dS e„ , (14) 

where e n is the unit normal direction at the surface. The integration is done over the whole surface. p m is the 
nuclear density and v is the average nucleonic speed obtained from the formula 

Hk 

v = V( — ^(Eav/a) 1 / 4 (15) 
nm 

with E av is the available energy and the level density parameter, a, is taken to be Acn/W- For the generalised 
shape (|la|), one-body friction, , is obtained as 

r+c 

7, OB = 2tt Pm vR 2 CN f(dc/dx) J dzp[l+ p/ 2 ]-^ 2 [A cP / + (1/2) pA'f (16) 

where pi, A' c are the derivatives of p, A c with respect to z and all other quantities are defined earlier. The tangential 
friction is calculated in a similar way as in Eqn. 

After the formation, the compound nucleus is in the minimum of the potential energy surface and it is assumed that 
the total initial energy available in the fusion process is equilibrated to allow the system to be described in terms of a 
thermodynamic state. Some part of it may be locked in rotational energy of the compound nucleus and the remaining 
part is the available excitation energy. However, what fraction of the available excitation energy will be converted 
into the collective excitation leading to the temporal evolution of fission shapes, is not known apriori. In our model, 
it is assumed that a random fraction of this available energy is imparted to the collective degrees of freedom which 
initiates the dynamics of the fission process. The initial conditions of r and r are 

r(t = 0) = r rn i n , (17a) 

r(t = 0) - {E^-Rn/Zp) 1 / 2 , (17b) 

where Rn is a random number between and 1 from uniform probability distribution. The available excitation 
energy, E*, is given by 

E* = E cm + Q eri t — L 2 /2Icn, (18) 

where E cm and Q en t are the centre of mass energy and Q-value in the entrance channel, and, Icn is the moment 
of inertia of the compound nucleus. The initial conditions of the angle variables at time t — are 

6{0) = 0i (0) - 9 2 (0) - 0'i (0) - 9 2 (0) = 0, (19a) 

0(0) = L/Icn (19b) 

The dynamical evolution starts at r = r m i n , where the minimum of the potential energy occurs. Once it reaches 
the saddle point or the top of the barrier, it is almost certain that it will reach the scission point. 
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C. The neutron multiplicity 



The emission of the prescission neutrons is incorporated in the present model as follows. During the temporal 
evolution of the fission trajectory the intrinsic excitation of the system is calculated at each time step. Correspondingly, 
the neutron decay width at that instant, r„, is calculated using the relation T n = HW n , where the decay rate W n is 
given by 

r E max ,2tt 

W -L dE iM- < 20 > 

The rate of decay A ^ A — 1 + ninan energy interval [E, E + dE] and a time interval [t, t + dt], iwf t , is given by 

dEdt = (1/?r H )Eam ^ UA (EX) ■ (21) 

The quantity loa(E*) is the level density for the nucleus with atomic number A and excitation energy E*. <Ji nv is 
the inverse cross section for the reaction (A — 1) + n — ► A and [i r is the reduced mass of (A—l,n) system. The upper 
limit of integration in Eqn. |2^ is given by 

E max = E* + Ba — [B A -i +B n ), (22) 

where Ba is the binding energy of the nucleus with atomic number A and B n is the neutron separation energy. 
The emission of neutron in the evolution of trajectory is now conceived by fixing a criterion. We evaluate the ratio of 
neutron decay time r n (= h/T n ) and the time step r of the calculation. The ratio r/r„ is compared with a random 
number Rn from a uniform probability distribution. The criterion of emission of a neutron at random time is fixed 
with the rule that whenever 



t/t„ > R N (23) 

the emission of a neutron takes place. If condition ( p3[ ) is not satisfied, no emission of neutron takes place. The 
probability of emission of a neutron in time r is (t/t„). The time step r is chosen in such a way that it satisfies the 
condition t/t„ <C 1. This suggests that the evaporation is a Poisson process leading to exponential decay law with 
half life r„ |22j| . Consequently the probability of emission of two or more neutrons in time r would be extremely small. 

The kinetic energy of the emitted neutron is extracted through random sampling technique. For this purpose, it is 
assumed that the system is in thermal equilibrium at each instant of time t, and therefore, the energy distribution of 
the emitted neutrons may be represented by a normalised Boltzmann distribution corresponding to the instantaneous 
temperature of the system. From a uniformly distributed random number sequence {x n } in the interval [0,1], we 
construct another random number sequence {y n } with probability distribution f(y), where f(y) ~ exp(— (3(t)y) is a 
normalised Boltzmann distribution corresponding to the temperature (3(t) at any instant of time t. Then, the sequence 
{y n } is obtained from the sequence {x n } by the relation, 

y(x)=F- 1 (x), (24a) 

where, F^ 1 is the inverse of the function F(y), which is given by, 

" f(y)dy = F(y). (24b) 

The inverse function is computed numerically by forming a table of integral values. The energy of the emitted 
neutron is given by E n — yE" Lax , where E" lax is chosen in such a way that the Boltzmann probability at that energy 
is negligible for all instants of time t. After the emission of the neutron, the intrinsic excitation energy is recalculated 
and the trajectory is continued. In this way for each angular momentum / of the compound system, the average 
number of emitted neutrons per fission event < M n is calculated. 




D. The statistics 



As argued previously, in the present model randomness is introduced only at the initial level when the compound 
nucleus is at the minimum of the potential energy surface. A random fraction of excitation energy is given to the 
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collective degrees of freedom. As a consequence, all trajectories would not be able to cross the barrier and would 
not have fission fate. Apart from that, there is a parameter a in Eqn. [la| which decides the final asymmetry of 
fission fragments uniquely through Eqn. 0. For obtaining different asymmetry, one should introduce a probability 
distribution Pi (a) at the initial level. It is natural to assume Pi {a) to be proportional to the density of states available 
for that a. Thus Pi(a) is taken to be 

Pi(a) oc u- 2 exp[2(au) 1/2 } (25a) 
where a is the level density parameter and u is given by 

u = E*-V min (a) (25b) 
and V m i n (pt) is the minimum of the potential energy surface for a given a. 

The compound nucleus is formed from the fusion process with different angular momentum. The dynamics to 
follow after its formation depends intricately on this angular momentum as seen from ([l8]). Hence the probability to 
cross the barrier would depend upon this angular momentum. We call this probability Pi(f, a\l). This is obtained as 
the ratio of number of trajectories crossing the barrier for given a and I and the total number of trajectories chosen. 
Therefore average of any observable quantity, O, is given by 

< Q >= E;;o r (2/ + l)0(a,OP ; (/,a|OF ; («) m 

E!=o er (2Z + lW,«|0fl(a) ' 

where the quantity O may be any of the rclcvcnt observablcs of interest, e.g., kinetic energy, neutron multiplicity 
etc., and, l cr is the critical angular momentum for fusion. 



III. RESULTS AND DISCUSSIONS. 



A large amount of prescission neutron multiplicity data over a wide range of excitation energies and masses of 
the compound nucleus is presently available in the literature. We have chosen a few representative systems in the 
ranges of masses A C n ~ 150-250, and, excitation energies Eq N ~ 60-160 MeV. All the systems considered here are 
above the Businaro-Gallone point and symmetric fission is the predominant mode of decay. Therefore, the theoretical 
predictions for various physical observables, ie, prescission neutron multiplicities (n pre ), total kinetic energies (TKE), 
average energies of the prescission neutrons (< E n >), etc., have been made for the the symmetric fission and 
confronted with the respective data. In addition, some recent experiments have been reported where fragment mass 
asymmetry dependence of the related physical observables have been studied for a few of the systems mentioned 
above. Explanation of such exclusive data is a crucial test for any theoretical model and it has not been, to the best of 
our knowledge, attempted so far. Therefore, we have made calculations for the fragment mass asymmetry dependence 
of some of the physical observables and the results have been compared with the corresponding experimental data. 



A. Fission shapes and Friction form factors 



The evolution of fission shapes for the symmetric (a = 0.0) as well as asymmetric (a = 0.15) splitting of a 
representative compound system of mass AcAf=200 are illustrated in Fig. [|. Fig. ||a to Fig. ||g represent the gradual 
evolution of the shape from the spherical ground state (c = 1) to the scission point (c = c sc ). It may be pointed out 
that the centre to centre separation at the scission point is ~ 18 fm which is much larger than what one gets in the 
schematic shape of p3| . 

In the present calculations, a combination of both one-body and two-body frictions has been used to calculate the 
fission trajectories and the fission observables. In Fig. ^ form factors of the one-body and two-body frictions are 
displayed as a function of the centre to centre separation between the two symmetric fragments for a typical system 
of mass Aqn — 200. It is seen from Fig. that at smaller separations ( when the shape is nearly mononuclear), 
one-body friction is stronger whereas at larger separations, two-body friction dominates. However, one-body friction 
does not change much with the increase in separation between the fragments. 
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B. Prescission Neutron Multiplicities 



As prescission neutron multiplicities depend on time scale of the fission process and vis-a-vis on the magnitude 
of nuclear friction, theoretical estimates of the friction coefficients are usually made by reproducing the prescission 
neutron multiplicity data using friction coefficient as an input parameter in the model. In the present calculations, one- 
body 'wall' friction has been used in the ground state to saddle region, where nuclear shapes are nearly mononuclear. 
The strength of the one-body friction used was attenuated to 10% of the original 'wall' value. This weakening of 
the wall friction has also been confirmed from the study of the role of chaos in dissipative nuclear dynamics p7[ . 
In the saddle to scission region, on the other hand, the nuclear dissipation was taken to be of two-body origin. For 



the two-body friction, the viscous drag was calculated in the framework of Werner- Wheeler using Eqn. (8a) and the 
value of the viscosity coefficient /io used in the present calculation was (4 x I0~ 23 Af eV ■ sec ■ /to -3 ). This value of /xo 
corresponds to 0.06 TP ( ITP = 6.24 x lO- 23 MeV ■ sec ■ fm' 3 ). 

The calculated prescission neutron multiplicities have been displayed in Fig. U as a function of the initial excitation 
energy of the compound nucleus for two different mass regions, ie, for Acn ~ 150 {upper half), and Acn ~ 200 (lower 
half). The solid curves represent the results of the present calculations and different symbols correspond to different 
sets of experimental data ||f|||. It is seen that for heavier systems (A C n ~ 200), the theoretical predictions are 
in good agreement with the corresponding experimental data. For lighter systems (Acn ~ 150), the experimental 
points are somewhat scattered and the theory is seen to reproduce quite well the average trend of the data. Here, the 
prescission neutron multiplicities are less and the uncertainties are more which are reflected in the larger error bars 
of the experimental measurements. Moreover, as different experimental points belong to different compound nuclei 
(different symbols in the figure), additional fluctuations in neutron emission due to specific structure effects may not 
be ruled out. For example, compound system 162 Yb, formed in the reaction 18 + 144 Sm (filled diamond), is quite 
neutron deficient compared to les Yb, formed in the reaction 18 + 150 Sm (open triangle). Therefore, neutron emission 
from the former is expected to be somewhat less. In fact, such system dependent fluctuations of the average neutron 
multiplicities have also been observed when calculations were done for specific systems (see text below) . At very high 
excitation energies, the observed multiplicity (open diamond) was found to be lower than the average theoretical trend. 
In this case, the incident energy was quite high (> 10 MeV/nucleon), and the onset of preequilibrium emission process 
may not be ruled out. As preequilibrium particles carry away a larger amount of energy compared to the evaporated 
particles, the fused composite cools down faster and subsequently leads to fewer emission of evaporated particles ||. 
As preeqilibrium emission has not been considered in the present calculation, the model predictions of multiplicities 
are expected to be somewhat higher than the experimental measurements of the same at higher bombarding energies. 

In Fig. ^(a), we have plotted prescission neutron multiplicities n pre as function of the compound nuclear mass 
Acn for 158.8 MeV ls O induced reactions on different targets. The solid curve represents the theoretical predictions 
whereas the filled circles correspond to the experimental data ||. It is seen that the observed multiplicities increase 
with the increase in Acn m general and show some fluctuations in the vicinity of Acn ~ 160— 170, in particular. The 
present calculations are found to be quite successful in reproducing the general trend of the data over the whole range 
of masses studied. The fluctuations in the observed multiplicities, which may be due to specific structure effects, as 
discussed earlier, have also been reproduced qualitatively in the present model. 

It is interesting to note that the present model with a single value of friction coefficient, [1q has been able to 
explain the general trend of the data over the whole range of masses and excitation energies of the compound nuclei 
studied. This is similar to what one expects in the Werner- Wheeler prescription, which predicts that the reduced 
friction coefficient, ft should be a universal function depending only on the collective degrees of freedom [^|. Earlier 
attempts, on the contrary, have shown that the values of the friction coefficients are system dependent and vary over a 
wide range (typically, the reduced friction coefficient, ft, may have values between 2 x 10 _21 sec _1 and 20 x I0~ 21 sec -1 
p3|). However, Frobrich et. al. attempted to arrive at such a universal value of the ft in the framework of modified 
Langevin equation model |C3 by using the level density formula proposed by Ignatyuk et. al. |2(| . They had to use 
a constant value of ft (= 2 X IO -21 sec -1 ) for the ground state to saddle followed by a set of ft values proportional to 
the elongation parameter, which reaches upto a value of 30 x IO -21 sec -1 at the scission point in order to reproduce 
the experimental data. 



C. Energy of emitted neutrons and TKE 

To have a closer look into the predictions of the present model so far as other related observables are concerned, 
the average energy of the prescission neutrons, < E n > and total kinetic energy (TKE) of the fragments have been 
plotted as a function of the compound nuclear mass, Acn, in Figs. ||(b) and ||(c), respectively. It is seen from the 
figure that the theoretical predictions for all the observables (solid curve) are in good agreement with the respective 
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experimental data (filled circles) . The observed system dependent fluctuations of both < E n > and TKE are also, 
reproduced fairly well in the present calculations. 



D. Fragment mass asymmetry dependence 

Calculations have also been performed to reproduce the available experimental data for the asymmetric fission 
of the heavy mass compound systems. The systems studied for the asymmetric fission are the ones produced in 
18 (Eiab — 158.8 MeV) induced reactions on 15i Sm, 197 Au and 238 U Figjfj] shows the predicted prescission 
neutron multiplicities n pre (dashed curves), calculated with the same friction values as used in symmetric fission 
in the previous section, as a function of the fragment asymmetry, alongwith the corresponding experimental data 
(solid circles). It is evident from the figure that, for all the systems, the predicted values of n pre are very weakly 
dependent on the fragment asymmetry, whereas the experimental values decrease rapidly with the increase in fragment 
asymmetry. It may be conjectured that the friction form factors calculated from the Werner- Wheeler prescription, 
which successfully explains the symmetric fission data, are rather inadequate in reproducing the experimental data in 
the case of asymmetric fission. It appears that the above friction form factors need some modification, ie, some extra 
asymmetry dependence should be included in the form factor to explain the asymmetric fission data. It has been 
found that, with the inclusion of a factor h(a) — exp(— Ka 2 ) in the expression for friction form factor (Eqn. |8a| ), the 
predicted values of n pre (solid curves) agree quite well with the respective experimental data. The value of the constant 
K was found to be 161 ± 3 which is independent of the mass of the compound system. It is, therefore, interesting to 
note that with the inclusion of this extra term h(a), we can still use the same value of the viscosity coefficient, /xq, 



as used earlier in Sees. IIIB, III C to explain the prescission neutron multiplicity data for both symmetric as well as 
asymmetric fission. 

We have also studied the variation of prescission neutron multiplicities as well as the time scale of fission as function 
of the angular momentum. The Figs. |?](a, b) show the calculated results for fission time Tf and neutron multiplicities 
npre, respectively, for a typical system ls O+ 197 Au for different values of asymmetry parameter a. It is clear from 
the figure that both n pre and Tf have only a weak dependence on L. However, both these quantities have a strong 
dependence on fragment asymmetry and are found to decrease with the increase in fragment asymmetry. The time 
scales vary in the range of 0.5 to 1.5 (xl0~ 20 sec.) as asymmetry parameter a decreases from 0.1 to 0. A similar 
trend is observed in the variation of n pre . 



IV. SUMMARY AND CONCLUSIONS 



To sum up, we have developed a dynamical model for fission where fission trajectories are generated by solving 
Euler-Lagrange equations of motion using a combination of one- and two-body frictions. Evolution of shapes of the 
fissioning nuclei have been computed from the generalised shape parametrisation of Brack et. al. ^ij. Different 
friction form factors have been used for ground state to saddle and saddle to scission regions. In the ground state to 
saddle region, the one body 'wall' friction has been used with an attenuation coefficient of 0.1 , whereas the two body 
dissipative forces, derived from Werner- Wheeler prescription, have been used in the saddle to scission region with the 
value of the two body viscosity coefficient u = 4 x lO~ 23 MeV ■ sec- fm~ 3 . Same values of the friction coefficients have 
been used for all the systems considered here. With these values of friction coefficients, the typical time scales of fission 
as obtained from the present calculations were ~ (1 — 2) x 10~ 20 sec. for symmetric fission, which are similar to the 
values reported earlier in the literature. Emission of neutrons along the fission trajectories has been simulated through 
Monte-Carlo simulation technique. The evolution of the fission trajectories has been corrected for prescission proton 
emission, which has been simulated in a similar way as it was done for neutron emission. However, the prescission 
proton emission was found to be quite small compared to the neutron emission for the systems considered in the 
present studies. The emission of complex particles ( d, t, 4 He etc.) and their effects on prescision neutron emission 
has been neglected in the present calculation as it is difficult to incorporate them in the present model because of their 
composite nature. Generally complex particle emission is much less as compared to the neutron emission; however, 
for some fissioning nuclei a emission may be relatively favoured due to structure effects and it may affect the neutron 
multiplicity and other fission observables. Theoretical predictions for prescission neutron multiplicities, total kinetic 
energies and the mean energies of the prescission neutrons, for all the systems considered here, agree quite well with 
the corresponding available experimental data. 

In the present model, the fragment mass asymmetry dependence of the prescission neutron multiplicities and other 
physically relevent quantities have also been studied. It has been found that with the inclusion of an extra factor 
(exp(— Ka 2 ), K—161±3) in the friction form factor (Eqn. pah, the observed decrease of the prescission neutron 
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multiplicities with the increase in fragment mass asymmetry, are explained fairly well with the same value of the 
viscosity coefficient fiQ as used in the case of symmetric fission. Further, it has been observed that the prescission 
neutron multiplicities as well as fission time are weakly dependent on the angular momentum of the fissioning system. 
However, the fission times are found to be strongly dependent upon the fragment mass asymmetry and decrease from 
1.5 to 0.5 (xl0~ 20 sec.) as the the asymmetry parameter a increase from 0.0 to 0.1. 

To conclude, in the present model, with the generalised shape parametrisation and modified Werner- Wheeler friction 
form factor, both symmetric and asymmetric splitting of the compound nucleus can be treated on the same footing. 
The present model is quite successful in explaining prescission neutron multiplicities, total fragment kinetic energies, 
average energies of the prescission neutrons, as well as their fragment mass asymmetry dependence. 
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FIG. 1. Variation of neck radius n with c for different values of the parameter a. Dash, solid and dash-dot lines correspond 
to a— 0, 0.05 and 0.10, respectively. 
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FIG. 2. Evolution of fission shapes for a compound system of mass Acn = 200 for a = 0.0 and 0.15. Figs, (a) to (g) 
represent the shapes of the fissioning shapes for c = 1, 1.2, 1.4, 1.6, 1.8, 2.0 and c sc , respectively. 'Negative' values of p 2 are 
the mirror reflection of the upper half about the elongation axis. 
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FIG. 3. Friction form factors for one-body and two-body frictions plotted as a function of the centre to centre separation r 
between the two symmetric fragments for a compound system of mass Acn = 200. 
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FIG. 4. Prescission neutron multiplicities plotted as a function of the initial excitation energy E* ni of the compound nuclei 
of masses Acn ~ 150 (upper half), and Acn ~ 200 (lower half). The solid curve is the present calculation. Different symbols 
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